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^ ; ABSTRACT 

We present high-resolution, iJ-band, imaging observations, conected with Subaru/HiCIAO, of the 
scattered hght from the transitional disk around SAO 206462 (HD 135344B). Although previous sub- 
, mm imagery suggested the existence of the dust-depleted cavity at r < 46 AU, our observations reveal 

the presence of scattered light components as close as 0'.'2 (~ 28 AU) from the star. Moreover, we 
have discovered two small-scale spiral structures lying within 0'.'5 (~ 70 AU). We present models for 
Ph ! the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h ~ 0.1, 

[t"! ■ which is consistent with previous sub-mm observations. This model can potentially give estimates of 

the temperature and rotation profiles of the disk based on dynamical processes, independently from 
sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable 
on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely 
Q ■ identify the origin of these spirals, planets embedded in the disk may be capable of exciting the 

| H I observed morphology. Assuming that this is the case, we can make predictions on the locations and, 

C/5 . possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths 

■ observations. 

Subject headings: circumstellar matter — instrumentation: high angular resolution — polarization 
, — protoplanetary disks — stars: individual (SAO 206462, HD 135344B) — waves 
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1. INTRODUCTION 

Dynamical processes in protoplanetary disks such as turbulence or disk-planet interaction are important in under- 
standing physical condition and evolution of disks, and planet formation processes. High resolution, direct imaging 
observations of circumstell ar / protoplanetary disk s can reveal non-axisymmetric structures, providing insight into such 
dynamical processes (e.g.. [Hashimoto et al.ll20lil ). 

Recent observations have identified a class of protoplanetary disks harboring tens of AU-scale holes/gaps at their 
centers: the so-called transitional disks. One well-studied system of that class is the rapidly rotating Herbig F star, 
SAO 20 6462 (HP 135344B, F4Ve, d = 142 ±27 pc, M = 1.7;°;iM q, IMiiUer et al.| [201 l|). The observation s of CO line 
profiles pent et al.ll2005HPontoppidan et aIll2008l:lLvo et al.l201ll ) and stellar rotation (jMiiller et al.ll2011[ ) consistently 
indicate an almost face-o n geometry (i ^ 1 1°)- The gap in the disk was predicted from the infrared (IR) spectral 
energy distrib ution (SEP . 'Brown et al.|[2007D. and wa s subsequently imaged in sub- mm dust continuum at ^ 0'.'5 x 0'.'25 
resolution (Br own et al., 2009). And rews et al.l (|201lD estimate the gap radius to be ^ 46 AU and the surface density 
within the gap to be 10~^'^ times smaller than that extr apolated from th e outer disk. The gas in Keplerian motion 
surrounding the gap regi on is also imaged by CO lines (iLvo et al.ll201lD. The C O rovibrational line observations 
(jPontoppidan et aLri2008| ) and [01] spectral line observations (jvan der Plas et al.l [20 081 indicate the presence of a 
gas disk in the vicinity (several AU-scale) of the star. SEP modeling (|Gradv et al.l i2009i) and NIR interf erometry 
(jFedele eT al. 2008) indicate the presence of an inner dust belt, which is temporally variable (jSitko et al.ll2012l ) and not 
coplanar with the outer disk (Benisty 2011, private communication). New imaging with high spatial resolution and 
sensitivity is required to understand the inner structures of the disk. The outer portions of gaps can n ow be resolved 
using 8-10 m ground-based telescopes at near infrared (NIR) wavelengths (e.g.. iThalmann et al.ll2010l for LkCa 15). 

In this Letter, we present iJ-band polarized intensity {PI) observations of the disk of SAO 206462 down to r ~ 
0'.'2 28 AU) scale at 0'.'06 (~ 8.4 AU) resolution. Interior to the sub-mm resolved gap, we find spiral structures, 
indicative of dynamical processes. We use the spiral density wave theory to interpret the structure, and estimate disk's 
physical parameters. 

2. OBSERVATIONS & DATA REDUCTION 
2.1. HiCIAO Observations 

SAO 206462 was observed in the iJ-band (1. 6 /im) using the high-contrast imaging instrument HiCIAO (iTamura et al.l 
120061: iHodapp et alll2008l: iSuzuki et ai]|2010[ ) on the Subaru Telescop e on 2 011 May 20 UT as part of Strategic Ex- 
plorations of Exop l anets and Pisks wi th Subaru (SEEPS. iTamural |2009[ ). The adaptive optics system (A0188; 
iHavano et all 120041: iMinowa et aLll2010[ ) provided a stable stellar point spread function (PSF, FWHM = 0'.'06). We 
used a combined angular differential imaging (API) and polarization differential imaging (PPI) mode with a field of 

view of 10 X 20 and a pixel scale of 9.5 mas pixeP^. A O'.'3-diameter circular occulting mask was used to suppress 
the bright stellar halo. The half-wave plates were placed to four angular positions from 0°, 45°, 22.5°, and 67.5° in 
sequence with one 30 sec exposure per wave plate position. The total integration time of the PI image was 780 sec 
after removing low quality images with large FWHMs by careful inspections of the stellar PSF. 

2.2. PDI Data Reduction 

The raw images were corrected using IRAFQ for dark current and fiat-field following the standard reduction scheme. 
We applied a distortion correction using globular cluster M5 data taken within a few days, using IRAF packages 
GEOMAP and GE OTRAN. Stokes {Q, U) parameters and the PI image were obtained in the standard approach 
(e.g., Hinklev et al.i.2009} as follows. By subtracting two images of extraordinary- and ordinary-rays at each wave 
plate position, we obtained +Q, — Q, +U, and —U images, from which 2Q and 2U images were made by another 
subtraction to eliminate remaining aberration. PI was then giv en by PI = y/ -\- U"^. Instrumental polarization of 
HiCIAO at the Nasmyth instrument was corrected by following iJoos et al.l ()2008[ ). From frame- by- frame deviations, 
the typical error of surface brightness (SB) was estimated to be ^ 5% at r 0'.'5 when averaged over 5x5 pixels (^ 
PSF scale). Comparing different data reduction methods (frame selections and instrumental polarization estimates), 
we expect that the systematic uncertainty of the SB of PI to be ~ 10%. 

2.3. Contemporaneous Photometry 

Since SAO 206462 shows variability in NIR wavelengths ()Sitko et al.l I20l3) . it is important to take photometry 
simultaneously with disk observations. 77-band photometry was obtained just before and after the disk imaging 
without the coronagraphic spot with the adaptive optics, by sixteen 1.5 sec exposures at four spatially dithered 
positions. An NPIO filter (9.8 ± 0.1% transmission) was used to avoid saturation. Using the MKO filter set, the 
iJ-band (Acff = 1.615 ^m, FWHM = 0.29 ^m, Tokunaga et al. 2002) magnitude was 6.96 ± 0.07 mag. 

Broad-band VRIJHK photometry was obtained on 2011 May 23-26, star ting within 48 hour s of the HiCIAO ob- 
servation, using the Rapid Eye Mount (REM) Telescope at La Silla, Chile (jCovino et al.ll2004[) . The REM i?-band 
filter has Aog = 1.65 /Ltm, FWHM = 0.35 /xm: broader and displaced to longer wavelengths than the MKO filter. The 
observed data were reduced differentially using SAO 206463 (AOV). The IR excess due to the inner disk (Figure [Ij was 

^ IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Fig. 1.— SED for SAO 206462 obtained by REM, adapted, in part, from lSitko et al.l l[20ll ). The REM observations consist of VRIJHK 
photometry. Also shown are spectra obtained with the SpeX spectrograph. See Section l2.3l for the data reduction techniques. 



average for the range observed in 2009-2011 ()Sitko et al.ll2012|) . No significant variation was observed during the 2011 
May observations, except for the smah long-term fading trend (Jm = .08 ± 0.02 mag) ove r the observation period. 

Figure [U also displays spectra obtained with the SpeX spectrograph (iRavner et al.ll200"9| ) on NASA's Infrared Tele- 
scope Facility (IRTF). The observations were obtained in the cross-dispersed (XD) ec helle mode between 0.8 and 5.1 nm. 
using a 0'.'8 slit {R - 900) and calibrated using HD 129685 (AOV) with SpeXtool (|Vacca et al.ll200l ICushing et alJ 
1200 41). The absolute flux calibration, to correct f or light loss at the spectrograph slits, was accomplished in two ways: 
using photometry and wide-slit spectroscopy (see lSitko et al.|[20T2l ). The March SpeX data were normalized using the 
REM photometry, obtained in the days immediately after the SpeX observations, and when the star was photometri- 
cally stable. In July, SAO 206462 and the calibration star were observed with the SpeX prism using a S'.'O slit, which, 
under good seeing and transparency conditions, produces absolute flux to 5% accuracy. The REM photometry at 
iJ-band in May is ^ 0.25 mag brighter than the Subaru data. The Subaru photometry shows a low value even when 
considering the filter difference, suggesting that the outer disk is illuminated efficiently. 

3. RESULTS 
3.1. Spiral Structure 

The SAO 206 462 disk can be tr aced in PI from 0'.'2 to - I'.'O (28 - 140 AU), similar to the range imaged with 
HST/NICMOS (|Gradv et al.ll2009[ ). but with a factor of - 4 greater angular resolution. The total PI is 9.87 mJy ± 
0.06% at 0'.'2 < r < I'.'O, which is 0.6% of the stellar intensity The total P I at 0'.'42 < r < I'.'O is 3.94 mJy ± 0.1% 
while the total intensity by HST/NICMOS FllOW is 9.7 mJy ()Gradv et al .|[2009,). The average SB of PI at r = 0'.'46 

is - 6 mJy/asec^, whereas the total intensity by HST/NICMOS F160W is 30 mJy/asec^ (|Gradv et al.ll2009D . Given 
the NICMOS data uncertainties, the polarization fraction is ^ 20 — 40%, assuming no PSF halo in the HiCIAO data 
and no variable self-shadowing/illumination in the disk. Our measured polarization fraction is similar to that of HD 

100546 fl4%+^gf°. [Quanz et al.l[20Tlh and AB Aur (- 25% - 45%, EerrineEaDllOOl). 

Figure [2] shows th e PI image. Th e region interior to 0'.'4 is not a void and we do not see clear structural evidence 
of the cavity wall in [Andrews et all (2011) model (i?cav — 46 AU ~ 0'.'33). We see spiral arcs SI (east) and S2 (south- 
west). The PI at the location of the spirals is ^ 30% larger than that extrapolated from the smooth outer profile 
(botto m of Figure 121) ■ The b rightest portions of the spirals r oughly coincide with the bright thermal emission peaks at 
12 nm (jMarihas et al.|[20Tll ) and lie inside the ring noted bv lDoucet et al.l ()2006l ). It is also noted that we see a dip in 
PI in the north-west, probably due to the depolarization in the minor axis direction (see below), and that we do not 
see large-scale, localized shadow that might be cast by the inner dust belt if highly inclined relative to the outer disk. 
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Fig. 2. — Top: PI image of SAO 206462 in the north-up configuration with log-stretch color scales. The filled orange circles at the center 
indicate the mask size (r = 0'.'15). The circles around have r = 0'.'2, exterior to which the features are considered to be real. The right 
panel is central region's closeup. Different color scales are used to enhance the spirals labeled as "SI" and "S2". The "Dip" may be due to 
depolarization. Bottom: PI profile along the red line in the top right panel. The arrow indicates the location of S2. The position errors 
are not shown for visibility. 

3.2. Azimuthal and Radial Profiles 

Here, we summarize the overall disk structure exterior to the spirals. Figure [3] (top panels) shows the azimuthal SB 
profiles at r = 0'.'5 and r = (X'T. At r ^ O'.'S, SB has maxima around position angle (PA) ~ 50° - 60° and 230° - 240°. 
Since the polarization is maximized at 90° scattering fe.g.. Graham et "all 120071 ), it is i mplied that the disk major 
axis is at PA ~ 50° — 60°, compa rable to estimates by CO observations: PA = 56° ± 2° by iPontoppidan et al.l ()2008f ) 
and 64° ± 2° bv lLvo et all (|20ll . We adopt PA = 55° for the major axis and z = 11° for inclination (see Section [1]). 
Our spiral model fitting results (the next section) are little aff'ected even if we assume a face-o n geometry. 

From CO observations, it is known that the south-west side is receding fe.g., lLvo et al.l[2011t ). Therefore, either the 
north-west or south-east side is the near side. We do not see an obvious forward scattering excess in the NIR image. 
However, since the spirals are typically trailing, it is inferred that the south-east is the near side. 

Figure [3] also shows the radial PI profiles along the major axis, wh ich is roughly consistent with in the outer 
part, indicating a flat (not highly flared) disk (jWhitnev and Hartmannl[l9 92). The radial slopes vary as PA from ^ —2 
to ~ -4.5 (fltting at 0'.'6 <r < I'.'O) or from ~ -2.5 to ~ -5 (fltting at O'.'S < r < O'.'Q), with shallower slopes typically 
appearing in the minor axis directions. However, r^^ is representative on average. This slope is observed in several 
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Fig. 3. — Top: Azimuthal SB profile at r = 0'.'5 (left) and r = 0'.'7 (right) with PA measured from north to east. Bottom: Radial profile 
along the major axis. Position errors indicate the FWHM of the PSF. SB errors are estimated using frame-by-frame deviations. 



other HAeBe disks (e.g., iFukagawa et al.ll2010l for total intensity data), although HD 97048 disk exhibits a shallower 
slope (|Quanz et al.ll2012D . 

4. SPIRAL STRUCTURE MODELING 

Among several features in the image, the most interesting one is the non-axisymmetric spirals. In order to understand 
them , we propose a model based on the spiral density wave theory (e.g., Lin and Shu 1964; Goldrcich and Trcmain^ 
[19781 119791 : fOgilvie and Lubowl 12002). assuming that NIR emission traces the disk surface density structure. With 
such a model, the spiral structures can be used to infer the disk temperature, independently of, for example, CO line 
observations. 

The shape of the spiral density wave is determined by the location of the launching point (corotation radius Tc) and 
disk's thermal and rotation profiles. When the disk rotation angular frequency is f2(r) cx r^" and the sound speed is 
c{r) oc , the shape of the wave far from rc is given by 



sgn(r - rc) 
70 X 



he 

1+/3 



1 + 13 l-a + (3 



(1) 



in the polar coordinate (r, 0), where — c{rc)/rcfl{rc) denotes the disk aspect ratio at r = and gives the phase. 
Equa tion ([1]) well approximates the shape of the density wave given by the WKB theory (jRafikov 2002; Muto et al.l 
1201 1[) . When the spiral is excited by a planet in a circular orbit, its location is '-^ (rc,6'o). Equation ^ has five 
parameters, {rc,9o,hc,a,f3). 

Two non-axisymmetric features, SI and S2, (Figure |4|) arc identified as follows. First, local maxima in the radial SB 
profiles normalized by (to take into account the dilution of the stellar fiux) are traced at every 1° step with data 
at 170° < PA < 360° (SI) and at 50° < PA < 190° (S2). The points near the minor axis (PA - 325°) are excluded 
because the structure there may be affected by depolarization. The points at PA > 200° (SI) and PA > 114° (S2) 
may be a part of axisymmetric rings since they are found to have constant radii. After removing these points, we have 
27 (SI) and 56 (S2) points as representing samples of non-axisymmetric spirals, with the opening angle of ^ 15° for 
both SI and S2. We estimate that the uncertainty of the location of the maxima is given by the FWHM of the PSF. 
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Fig. 4. — Left: Red and green points indicate the locations of the maxima of radial profiles for SI (top) and S2 (bottom). The green ones 
are used for the fitting by Equation |[T} (Section |4]|. Dashed lines indicate circles with r = 0'.'39 (top) and r = O'.'S (bottom). Right: The 
spiral shape given by Equation {ij with the best-fit parameters of SI (top) and S2 (bottom). Crosses show the locations of {rc,0o)- 

In order to fit t he non-axi s ymme tric structures by Equation ([T]), we fix a and /3 at 1.5 (Kepler rotation) and 0.4, 
respectiyely, as in iLyo et al.l ( 2011h . while other parameters are varied as (O'.'l < Tc < 0'.'9,0 < 9q < 27r,0.05 < he < 
0.25). Note that different yalues of (3 yield similar results. Since it is difficult to fit SI and S2 simultaneously, they are 
fitted independently. 

The "best-fit" parameters are (rc,6'o,/ic) = (0'.'39, 204°, 0.08) for SI (reduced = 0.52) and (rc,6'o,/ic) 
(0'.'9, 353°, 0.24) for S2 (reduced = 0.31). The spiral shapes with these parameters are shown in FigurelH However, 
the parameter degeneracy is significant. Figure [5] shows the parameter space of (rc^o) with 63.8% confidence level 
for he — 0.1 and he = 0.2. Note that in Figure[5l the "best-fit" of (rc, 6*0) is outside the domain of confidence in some 
cases because he is not the same as the best-fit. Despite the parameter degeneracy, the values of the aspect ratio 
which fit the shape of the spiral {he ~ 0.1) are consistent with those obtained from the sub-mm map of the disk (e.g., 
h ^ 0.096(r/100AU)°-i^: [A^rews et al.l[20Tl . 

The spiral density wave theory predicts that the pattern speed deviates from the local Kepler speed; 
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Fig. 5. — Parameter degeneracy of the fitting for SI (top) and S2 (bottom) with he = 0.1 (left) and he = 0.2 (right). The red points 
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{'>'c,6o) for the best-fit parameters. Note that the best-fit values of he are different from those shown in the figure. 

is not necessarily equal to ^{r). When rc — 0'.'5(~ 70AU), the spiral will move ~ 10° in a decade, corresponding to 
movement of O'.'l. Considering the PSF scale of our observations and the locations of the spirals, such deviations can be 
detectable over a couple of decades. Moreover, if the two spirals have distinct corotation radii, their relative locations 
change in time due to the pattern speed difference. Such measurements will confirm that the observed feature is really 
the density wave, providing indisputable evidence of dynamical activity. 

Note that it would be difficult to detect spirals in colder disks (smaller h^), where spirals are more tightly wound, 
due to the blurring by the PSF. The lower detectable limit of he is typically he ^ 0.01 — 0.03 for our set of parameters. 
The combination of high angular resolution and warm temperatures allows the spiral structure in the SAO 206462 
disk to be resolved. Further spirals might be detectable in similarly warm disks. 

5. SUMMARY AND DISCUSSION 

In this Letter, we present a high-resolution image of the SAO 206462 transitional disk using Subaru/HiCIAO, with 
an inner working angle of 0'.'2. We discover non-axisymmetric spiral features, which can be explained by the spiral 
density wave theory with a reasonable value of the disk aspect ratio {he ^ 0.1). The model is robust in a sense that 
it does not assume the origin of such structure explicitly. 

The detection of scattered light within the sub-mm cavity itself is interesting, since [Andrews et al.l (|2011[ ) predicts 
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that the sub-mm cavity is heavily depleted. Our d ata in tandem wit h the millimeter data may suggest that the 
depletion of grains at different sizes is not uniform. iDong et all (|2012D discuss such discrepancies between sub-mm 
dust continuum images and NIR scattered light images in terms of differing spatial distributions as a function of grain 
size from a general theoretical perspective. 

Our major assumption is that PI at _ff-band, tracing the scattered light at the disk surface, actually traces the surface 
density variations. This assumption is valid when the d isk is in vertical, isothermal hydrostatic equilibrium without 
rapid radial surface density variations fe.g.. [Mutcill20lil ). Structures near the midplane are, however, preferentially 
observed at longer wavelengths with high spatial resolution: Atacama Large Millimeter /Submillimeter Array (ALMA) 
can be an ideal instrument. 

Among several possible causes for the spiral structures (see also iHashimoto et al.ll20l"ll ). one interesting idea is that 
planets excite them. In this case, the domain of possible locations of the planets is given by the green curves in Figure 
[5l If the two spirals have distinct corotation radii, there may be two (unseen) planets embedded in the disk. The 
amplitude of the surface density perturbation scales with the planet mass as (5S/S ~ GMp^/c^ for non-gap-opening 
low-mass planets (GMpVl/c^ <^ 1, e.g., iTanaka et al.ll200^ . In our data, the amplitude of the spiral wave is typically 
- 30% (Figure E]), implying Mp ~ 0.5Mj. 

The typical error of PI in our image is ^ 5%. If this is typical of HiCIAO, it is capable of detecting the indirect 
signatures of planets down to Mp ^ 0.05Afj. ADI is promising in finding a point source; however, small field rotation 
(^ 13°) due to the southerly declination (—39°) of SAO 206462 makes obtaining the total intensity difficult in our 
data. Future L-band observations may reveal thermal emission from a planet, if it exists, or its surrounding (accreting) 
gas. 
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18540238 (SI), and 22-2942 (TM), WPI Initiative, MEXT, Japan (ELT), NSF AST 1008440 (CAG), 1009203 (JC), 
and 1009314 (JPW), and NASA NNH06CC28C (MLS) and NNX09AC73G (CAG and MLS). Part of this research was 
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